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Abstract. A general, uniform, rigorous and constructive thermodynamic ap- 
proach to weakly nonlocal non-equilibrium thermodynamics is reviewed. A 
method is given to construct and restrict the evolution equations of physi- 
cal theories according to the Second Law of thermodynamics and considering 
weakly nonlocal constitutive state spaces. The evolution equations of internal 
variables, the classical irreversible thermodynamics and Korteweg fluids are 
treated. 



1. Introduction 

Weakly nonlocal, coarse grained, phase field and gradient are attributes of the- 
ories from different fields of physics indicating that in contradistinction to the tra- 
ditional treatments, the governing equations of the theory depend on higher order 
space derivatives of the state variables. The origin of idea goes back to the square 
gradient model of van der Waals for phase interfaces pQ, where it is extensively 
applied 13 HI [S] . Later applications go far beyond phase boundaries or thermo- 
dynamics. Nowadays weakly nonlocal is a nomination in continuum physics dealing 
with internal structures [3 El [10l Qj] , coarse grained or phase field appears in 
statistically motivated thermodynamics [13 UHl HH H21 [5] , and gradient is frequently 
used in mechanics in different context [Ta[T7l[ia[TS|5Dlinil21ll3IMI23liai23. 
The simplest way to demonstrate the meaning of weakly nonlocal extensions can be 
exemplified by the Ginzburg-Landau equation which is not only a specific equation 
in superconductivity, as it was introduced originally by Landau and Khalatnikov 
[28j , but a first weakly nonlocal extension of a homogeneous relaxation equation of 
an internal variable. The traditional derivation of the Ginzburg-Landau equation is 
based on a characteristic mixing of variational and thermodynamic considerations. 
One applies a variational principle for the static part and the functional derivatives 
are introduced as thermodynamic forces into a relaxation type equation. A clear 
variational derivation to obtain a first order differential equation is impossible with- 
out any further ado (e.g. without introducing new variables to avoid the first order 
time derivative, which is not a symmetric operator) [29) . One can apply these kinds 
of arguments in continuum theories in general, preserving the doubled theoretical 
framework separating reversible and irreversible parts of the equations E01 EH • 
However, there are also other attempts to unify the two parts with different addi- 
tional hypotheses and to eliminate this inconsistency of the traditional approach 

G21G2liB]. 

The ultimate aim is to find a unified, general, rigorous and predictive theoretical 
framework that makes possible to extend the governing equations of physics with 
higher order gradients of the continuum fields, beyond the traditional terms. The 
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method should be uniformly applicable from classical systems in local equilibrium 
up to relativistic systems beyond local equilibrium; should be general to incorpo- 
rate most of the mentioned classical examples of weakly nonlocal theories without 
specific assumptions; should reduce the independent additional assumption to a 
minimum and should be constructive to give calculational methods for systematic 
higher order extensions of the constitutive space. 

This paper is a general tutorial to the mathematical framework of such a theory. 
In the first section a general methodology of exploiting the Second Law for weakly 
nonlocal systems is given. The Second Law is considered as a constrained inequality 
where the constraints are the evolution equations of the system and their deriva- 
tives, depending on the order of the nonlocality. At the third section evolution 
equations of internal variables and their different weakly nonlocal extensions are 
treated. A first order weakly nonlocal theory leads to relaxation type ordinary dif- 
ferential equations, a second order nonlocality leads to Ginzburg-Landau equation 
and a second order nonlocal theory to dual internal variables unifies the evolution 
equations of internal variables derived by mechanical methods (by variational prin- 
ciples and dissipation potentials) and by thermodynamics (by heuristic application 
of the Second Law). At the fourth section we show that classical irreversible ther- 
modynamics can be incorporated naturally in our treatment, a first order weakly 
nonlocal theory of balance type evolution equations leads to the thermodynamic 
flux-force relations of classical irreversible thermodynamics with gradients of the 
intensives as thermodynamic forces. Finally we demonstrate the applicability of 
the method to one component heat conducting Korteweg fluids, that are first order 
weakly nonlocal in the energy and in the velocity and second order weakly nonlocal 
in the density. In that case nontrivial forms of the pressure tensor ensure the com- 
patibility to the Second Law. As a particular example we derive the constitutive 
functions of the Schrodinger-Madelung fluids. Finally a summary and discussions 
follow. 

2. Second law and weakly nonlocal constitutive spaces 

In this section we shortly summarize some methodological specialities of exploit- 
ing the Second Law in weakly nonlocal systems. One can find some further details 
in [35l [36] , where the role of the entropy flux is treated and in [37] , where the key 
idea is introduced, the derivative of the constraints as additional constraint. 

There are several interpretations of the Second Law in non-equilibrium thermo- 
dynamics (see e.g. [35] [3^]). If we wanted to ensure that the Second Law will be a 
consequence of material properties, then we should postulate that the entropy was 
an increasing function along the solutions of the evolution equations. With this 
assumption we are compatible with the classical heuristic method that puts the 
evolution equations in the balance of the entropy, and exploits the consequences 
of the inequality. In the classical approach a quadratic expression is recognized 
and then a relation is sought between the introduced thermodynamic fluxes and 
forces [IHl HI] • Here the evolution equations are considered as constraints for the 
inequality of the entropy production [42]. We arrive to the Coleman- Noll proce- 
dure recognizing the algebraic part of the problem where the different derivatives 
are independent. When instead of putting constraints into the evolution equations, 
the constrained algebraic inequality is solved by multipliers then the calculation 
is called Liu procedure |43j . The Coleman-Noll and Liu procedures are equivalent 
in simple systems [44], but the later one preserves the symmetries of the evolution 
equations and the constraints. Liu procedure is based on a linear algebraic theorem, 
called Liu's theorem in the thermodynamic literature [451 46j and an interpreta- 
tion of the role of entropy inequality. Hauser and Kirchner recognized that Liu's 
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theorem is a consequence of a famous statement of optimization theory and linear 
programming, the so called Farkas's lemma [47] . That theorem was proved first by 
Farkas in 1894 and independently by Minkowski in 1896 [48]. In the Appendix we 
have formulated and proved the classical Farkas lemma, the affine Farkas lemma 
and Liu's theorem. 

Originally Farkas developed his lemma to formulate correctly the Fourier's prin- 
ciple of mechanics of mass-points, which is the generalization of d'Alembert's prin- 
ciple in case of inequality constraints [49 . The role of Liu's theorem in continuum 
physics is similar in some sense: we want to give the correct form of the evolution 
equations taking into account the requirement of the Second Law, that is the en- 
tropy inequality. In fact only the material part of the evolution equations - the 
constitutive functions - is/are restricted. The mathematical formulation introduces 
a kind of dual point of view, because in Liu procedure the evolution equations - the 
partial or ordinary differential equations determining the evolution of the system - 
is a condition, a constraint for the entropy inequality. The variables (e.g. fields) 
in the evolution equation form the basic state space. The constitutive quantities 
depend on these functions, on the basic state and some of its derivatives. To make 
the problem algebraically manageable the basic state variables and their deriva- 
tives are considered as independent quantities. Some of them can be incorporated 
into the constitutive state space (or large state space [46]), into the domain of the 
constitutive functions. The entropy inequality, our objective function, has a special 
balance form, and determines the process directions, the independent variables of 
the algebraic problem: those are the derivatives of the constitutive state, that are 
not included in the constitutive space. The choice of the constitutive state space is 
crucial and can result in different kind of restricted constitutive functions with Liu 
procedure. 

As we have already mentioned, the weakly nonlocal constitutive spaces contain 
space derivatives. In this cases some derivatives of the constraints - balances and 
others - can appear as additional constraints that further restrict the entropy in- 
equality. Whether these additional constraints should be considered or not depends 
on the order of weak nonlocality, the structure of the constraints and physical con- 
siderations. This is the peculiarity of the exploitation of the Second Law for weakly 
nonlocal systems. 

In the following sections we will give several examples to demonstrate the appli- 
cation of the formalism. 

3. Thermodynamic evolution of internal variables 

In this section we investigate the thermodynamic restrictions on the evolution 
equation of internal variables in continua at rest. First in a first order weakly 
nonlocal constitutive state space, then in a second order weakly nonlocal one and 
finally considering the peculiarities of dual internal variables. Further details of the 
related calculations and the physical interpretation can be found in [37] [501 EI] • 

3.1. First order nonlocality - relaxation. Let us investigate the thermody- 
namic restrictions for the evolution of a classical internal variable field a(t, r), in 
a continuum at rest related to an inertial observer. There are no constraints or 
knowledge regarding the form of its evolution equation. Therefore the evolution 
equation can be given in a general form as 

(1) dta + f = 0, 

with an arbitrary constitutive function /. Here and in the following the partial time 
derivatives are denoted by dt and the constitutive quantities are denoted by a hat 
("). First we introduce a first order weakly nonlocal constitutive state space spanned 
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by the basic fields a and their gradients c^a, where i = 1,2,3. The notation with 
indices with Einstein summation convention is applied by distinguishing covariant 
and contravariant components (vectors and covectors) by upper and lower indices, 
e.g. diJ l {— DivJ — V • J) denotes the divergence of the vector field J\ 

The above evolution equation is not completely arbitrary, it is restricted by the 
Second Law of thermodynamics. Therefore we assume that there is an entropy 
balance with a nonnegative production term 

(2) d t s + d l J l >0. 

Here the entropy density s and the entropy flux J 1 are constitutive quantities, 
too. Therefore in this case 

- the basic state space is spanned by a, 

- the constitutive state space is spanned by (a, did), 

- the constitutive functions are s, J 1 and /. 

We may develop the derivatives according to the constitutive assumptions as 

(3) d a sd t a + dg ia s dua + d a J l did + d dja J l d^a > 0. 

Here the partial derivatives of the constitutive functions are denoted in an ab- 
breviated manner, e.g. g| = d a s. The underlined partial derivatives of the basic 
field are not in the constitutive space, they are independent algebraic quantities 
and span the process direction space (<9(a, d^a, c^a). Then we can apply the the- 
orem of Liu identifying the underlined partial derivatives in ([1]) and © by p and 
their coefficients by a and b respectively. Now we introduce a Lagrange-Farkas 
multiplier A, which is a constitutive quantity, for the evolution equation JT]) and 
apply Liu procedure to determine the form of the constitutive functions, required 
by the entropy inequality: 

(4) < (d a s - Xj dta + d dza § dua + d dja J l djjCi + d a J l - A/. 

Here the existence of the multiplier A follows from Liu's theorem and we will 
exploit that the process directions are independent of the constitutive space in 
the sense that the values of these derivatives can be different for same values of its 
coefficients depending e.g. on the initial conditions of (JTJ) . Therefore the multipliers 
of the underlined terms in © are zero and give the Liu equations: 

(5) d t a : d a s = A, 

(6) d it a : d dia s = Q l , 

(7) d ija : <)„„.!■ 0". 

The first equation determines the Lagrange-Farkas multiplier and the last two 
ones show that both the entropy and the entropy flux is local, independent of the 
gradient of a. This is a complete solution of the system. The residual inequality is 

(8) 0<dj l (a)~d a s(a) f(a,d ia ). 

Now we assume that the residual entropy flux is zero J 1 = 0*. This is the 
situation in isotropic materials according to the representation theorem of isotropic 
vector functions. Therefore the entropy inequality reduces to a flux-force system. 
If we want to determine the form of the evolution equation then the entropy should 
be considered as a given function and / is an undetermined constitutive quantity. 
The classical solution of the above inequality is that the constitutive function / is 
proportional to the entropy derivative by a nonnegative constitutive multiplier: 

(9) f = -ld a s, f>0. 
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This is a general solution of the inequality if we assume two times differentiability 
of / (due to the mean value theorem of Lagrange, see e.g. [33]). In the practice 
we rarely exploit this generality and restrict ourselves to the constant I case. This 
remark applies further in ([15)). dSSJ)- (|37J) and ([55)) . 

Therefore the evolution equation of an internal variable restricted by the Second 
Law in isotropic materials is a relaxation type differential equation: 

(10) d t a = id a s. 

Now let us summarize the most important steps of the procedure: 

- Identification of the basic inequality and the basic constraints and the do- 
main of the corresponding functions (constitutive state space). 

- Performing the partial derivations and the identification of the process di- 
rection space by the partial derivatives of the basic state space and intro- 
ducing the additional (derivative) constraints in case of necessity. 

- Application of Liu's theorem. 

- Solution of the Liu equations and the dissipation inequality. 

In this example there was no need to apply additional derivative constraint, but 
it can be necessary incase of more extended constitutive state spaces as one can see 
in the next problem 

3.2. Second order nonlocality - Ginzburg-Landau equation. In this case we 
face to a similar problem, we are to determine the thermodynamic restrictions on 
the general evolution equation (fT]) , but now we assume that there is a second order 
weakly nonlocal state space spanned by the basic field a and its first and second 
space derivatives did and 9ya. 
Therefore in our second example 

- the basic state space is spanned by a, 

- the constitutive state space is spanned by (a, d%a, di a), 

- the constitutive functions are s, J 1 and /. 

The corresponding process direction space is spanned by the next derivatives 
(dta,dua,duja,dijka). Let us observe that these are not independent any more, 
the gradient of |T]) is a linear relation on the process direction space. Therefore we 
should consider 

(11) d u a + d t f = t 

as a further constraint to the entropy inequality (0). We introduce the Lagrange- 
Farkas multipliers A for (fT)) and A* for (TTT)) . Let us apply Liu procedure again, but 
in this case not separating the different steps: 

< dts + dj 1 -\(d t a + -^(dua + dif) 

= d a s d t a + dg z aS d tt a + dg zja s d iot a + d a J l d t a + dg ]a J % %a + do ]ka J l d i3k a 

- A (d t a + /) - A 4 {dua + d a f did + dd 3 af dijd + dd Jk af dijkoj 
= (d a s - A J dta+ (d dia s - AM d^a + d dtja s d ijt a + 

{d d]h aJ l - k l d 9]ka f) d t3k a + d a J l d t a + d dja J l %a - 

(12) A 1 (d a f d t a + d dj J %a) - A/. 

We can see that the degenerations are different than in the previous case. In the 
following we do not give the detailed expositions of the theorem, but everybody 
can reconstruct that from the detailed presentation of the calculations. 
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The multipliers of the underlined partial derivatives in (fT2"|) . the process direction 
space, give the Liu equations: 



(13) d t a : d a § = A, 

(14) d lt a : d di J = A i , 

(15) d ljt a : d 9ija s = ij , 

(16) d t3k a : d dika J l = A l d d]ka f. 



The first two equations determine the Lagrange-Farkas multipliers by the en- 
tropy derivatives and the third equation shows that the entropy is independent of 
the second gradient of a. As a consequence, the Lagrange-Farkas multiplier A 4 is 
independent of that derivative, therefore the last equation can be integrated and 
gives 

(17) J l (a,dia,dija) = do ia s(a,dia) f(a, did, dijo) +T(a, did), 

where the residual entropy flux, 3\ is an arbitrary constitutive function and the 
variables of the constitutive functions are explicitly written. This is a complete 
solution of the system of the Liu equations ([T3)) - (fTo| . Therefore the dissipation 
inequality simplifies to the following form 

(18) < ditf + {di{d dia s) - 8 a s) /. 

Assuming that the residual entropy flux is zero y = 0\ the entropy inequality 
reduces to a flux-force system. In this case the classical solution of the inequality 
is 

(19) f = l(di{d 9ia s)-d a s), l>0. 

Therefore the form of the evolution equation of an internal variable in a second 
order weakly nonlocal constitutive state space is the Ginzburg-Landau equation: 

(20) d t a = l (d a s - d t (d dta s)) . 

We can get the classical form of the equation with a particular form of the entropy 
density s(a,dia) = s a (a) — ~f(dia) 2 , with quadratic dependence on the gradient and 
with a constant coefficient 7. 7 > because of the concavity of the entropy. Other 
nonlocal thermodynamic potentials also can be used in the above derivation, e.g. 
the free energy, but one should be careful to use a correct thermodynamic structure 
when coupling to other thermodynamic interactions beyond the one related to the 
internal variable [13] . 

Ginzburg-Landau equation and its variants appear in different fields of physics 
and are applied to several phenomena. Beyond their original appearance in super- 
conductivity they play an important role in pattern formation and they are the 
prototypical phase field models [HI [TS] . 

As we have mentioned in the introduction the traditional derivation of the 
Ginzburg-Landau equation has two main ingredients: 

- The static, equilibrium part is derived from a variational principle. 

- The dynamic part is added by stability arguments (relaxational form). 
The physical content of the two ingredients of the classical derivation is sound 

and transparent [T^]. On the other hand, the origin of the variational principle, 
the coupling of the two parts and the role of the Second Law of thermodynamics is 
ad-hoc and is not compatible with the general balance and constitutive structure of 
continuum physics. Here we unified the ingredients in a thermodynamic derivation, 
and we derived the form of the entropy flux, too. We did not refer to any kind 
of variational principle, however, the derived static part has a complete Euler- 
Lagrange form. The dynamic part contains a first order time derivative, therefore 
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one cannot hope to derive it from a variational principle of Hamiltonian type [23] . 
The static part belongs to the zero entropy production, in this sense it is reversible. 
In our approach we get the "reversible" part as a specific case of the thermodynamic, 
irreversible thinking. 

There are several alternate derivations based on different concepts [33J IE] • Here 
we have demonstrated, that the physical assumptions to get the Ginzburg-Landau 
equation are very moderate, the typical Ginzburg-Landau structure is a straight- 
forward consequence of the entropy inequality without any further additional set 
of concepts like the microforce balance or a variational principle. 

3.3. Dual internal variables - Hamiltonian structure. There are two classical 
approaches to determine the evolution of internal variables. 

When the evolution equations of internal variables are constructed exploiting 
the entropy inequality, using exclusively thermodynamic principles, then the corre- 
sponding variables are called internal variables of state [5]. This frame has the 
advantage of operating with familiar thermodynamic concepts (thermodynamic 
force, entropy), however, no inertial effects are considered. The thermodynamic 
theory of internal variables has a rich history (see the historical notes in [55]). A 
first more or less complete thermodynamic theory was suggested by Coleman and 
Gurtin [53] , and the clear presentation of the general ideas of the theory was given 
by Muschik [53] . Internal variables of state were applied for several phenomena in 
different areas of physics, biology, and material sciences. A complete description 
of the thermodynamic theory with plenty of applications based on this concept of 
internal variables of state can be found in [55] . 

There is a second method that generates the kinetic relations through the Hamil- 
tonian variational principle and suggests that inertial effects are unavoidable. This 
approach has a mechanical flavor, and the corresponding variables are called dy- 
namic degrees of freedom. Dissipation is added by dissipation potentials. This 
theoretical frame has the advantage of operating with familiar mechanical concepts 
(force, energy). The method was suggested by Maugin 56J, and it also has a large 
number of applications [57] [55] . The clear distinction between these two methods 
with a number of application areas is given by Maugin and Muschik [55] |5j|] and 
Maugin [9]. 

Here we follow the terminology of Maugin and Muschik [55] with some important 
extensions. We call internal variables of state those physical field quantities - beyond 
the classical ones - whose evolution is determined by thermodynamical principles. 
We call internal degrees of freedom those physical quantities - beyond the classical 
ones - whose dynamics is determined by mechanical principles. 

One of the questions concerning this doubled theoretical frame is related to 
common application of variational principles and thermodynamics. Basic physical 
equations of thermodynamical origin do not have variational formulations, at least 
without any further ado [2§]. That is well reflected by the appearance of dissipa- 
tion potentials as separate theoretical entities in variational models dealing with 
dissipation. 

On the other hand, with pure thermodynamical methods - in the internal vari- 
ables approach - inertial effects are not considered. Therefore, the coupling to sim- 
plest mechanical processes seemingly requires some additional assumptions, those 
are usually new principles of mechanical origin. 

On the other hand, with pure thermodynamic methods in the internal variables 
approach inertial effects arc not considered. Therefore, the coupling to simplest 
mechanical processes seemingly requires introducing some improvements, which are 
usually new principles of mechanical origin. 
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In the following we show that the mechanical structure arises from thermody- 
namic principles. Our suggestion requires dual internal variables and a particular 
generalization of the usual postulates of non-equilibrium thermodynamics: we do 
not require reciprocity relations. With dual internal variables we are able to get 
inertial effects and to reproduce the evolution of dynamic degrees of freedom. It 
could be impossible with a single internal variable. This is the price we pay for the 
generalization. In other words, instead of the doubling of the theoretical structure 
we suggest the doubling of the number of internal variables. 

Let us consider a thermodynamic system where the state space is spanned by two 
scalar internal variables a and b. Then the evolution of these variables is determined 
by the following differential equations 

(21) dta + f = 0, 

(22) d t b + g = 0. 

The functions / and g are constitutive functions, restricted by the Second Law 
of thermodynamics. The entropy inequality, the main ingredient of the Second Law 
can be written in the same form as previously in ([2]). The domain of the constitutive 
functions (our constitutive space) is spanned by the state space variables and by 
their first and second gradients. Therefore in this case 

- the basic state space is spanned by (a, 6), 

- the constitutive state space is spanned by (a, did, dija, b, dib, dijb), 

- the constitutive functions are s, J 1 , / and g. 

This is a weakly nonlocal constitutive space with second order weak nonlocality in 
both variables. The corresponding process direction space is spanned by the next 
derivatives {dta, dua, du 3 a, dij k a, dtb, dub, dujb, dijkb). The gradients of the evolu- 
tion equations ([2"Tj) - (j2"2"|l are constraints for the entropy inequality in the framework 
of a second order constitutive state space for both of our variables 

(23) 8 u a + dif = 0i, 

(24) dub + d ig = 0i. 

We introduce the Lagrange- Far kas multipliers A a , X b for ([2"Tj)-([2"2l and A* , A^ for 
(|23l - (|24|) . respectively. Liu procedure results in the Ginzburg-Landau structure of 
the previous subsection in a doubled form 

< dts + diJ'-Xafdta + fy -Ai(d ti a + dJ^ - 

\ b (d t b + g)~Al {d ti b + dig) 
= d a s d t a+ da z aS d it a+ dg ija s d ijt a+ d a J l d t a+ do ja J l d i3 a+ da ]ka J l d^a- 

A a (d t a + /) - A* (d ti a + d a f did + d 9j af d i3 a + d d]ka f d l]k a+ 
dbf dib + d 9j bf dijb + d 9jk bf d ijk bj 
+ d b s d t b + d 8lb s d it b + d dljb s d ijt b + d b J l d t b + dg ]b J l d^b + dg ]kb J l d l]k b - 
Ah {d t b + g)- k\ {d tl b + d b g d t b + d djb g %6 + d djkb g d ijk b+ 

(25) d a g did + d dja g d tJ a + d djka g d ijk a) . 



After some rearrangements one can get 
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< {d a s - X a )d t a + (d dia s - k l a )d lt a + d dijC1 s d ljt a + 
{d d]k aJ l - Kdd ]k af - kid djha g)d ijk a - 
XJ - (dj d t a + d 9]a f d ija + d b f d t b + dg jb f d i:j bj + 

d a J l dia + d d]a J l d lj a 
+ (d b s - X b )d t b + (d 9ib s - kl)d it b + d di]b s d i]t b + 
(do^bJ 1 - K9d jk bf - Kd 9jkb g )d ljk b - 
X b g - k\ (d b g dib + d djb g d i:j b + d a g d t a + d dja g d tj a) + 
d b J l dib + d djb j l dijb. 
Here the multipliers of the process direction space give the Liu equations: 



(26) 


d t a : 


: d a s — A Q , 




(27) 


d it a : 


: d dia s = Ai, 




(28) 


d t b : 


: d b s = A h , 




(29) 


d it b : 


: d dlb s = Al, 




(30) 


d ljt a : 


: <h, „s II''. 




(31) 


dijtb : 


: <)r, ()'•'. 




(32) 


d ljk a : 


: d djka f = Kd djka f- 


f Hd 9 . ka g 


(33) 


dtjkb : 


■ ddjubJ* = H,9a jk b9 4 





The first four equations determine the Lagrange-Farkas multipliers by the en- 
tropy derivatives. The fifth and the sixth one show that the entropy is independent 
of the second gradient of a and b. Consequently, the Lagrange-Farkas multipliers 
A l a and A^ are independent of that derivatives, therefore the last two equations can 
be integrated and give 

(34) J i = d dza s f + d dlb s g + y(a, d t a, b, d t b). 

Here the variables of the the residual entropy flux J 1 , an arbitrary constitutive 
function, are explicitly written. This is a complete solution of the system of Liu 
equations (f2"B" ]) - (|33p . Therefore the dissipation inequality simplifies to the following 
form 

(35) < da' + (di(d dia s) - d a s) f + {di{d dib s) - d b s) g. 

Now assuming that the residual entropy flux is zero, Z l = 0', the entropy in- 
equality reduces to a flux-force system of the following form: 

a-force: A = di{do ia s) — d a s, a-flux: /, 

b-force: B = di(do ib s) — d b s, b-flux: g. 

The classical solution of the entropy inequality gives a coupling of fluxes and 
forces as a system of coupled Ginzburg-Landau equations: 

(36) d t a = f = h (dtidd.aS) - d a s) + hi (di(d dib s) - d b s) , 

(37) d t a = g = l 2 i (di(d 9ia s) - d a s) + h [di[d d . b s) - d b s) . 

The constitutive Onsagerian coefficients h,h, li2, hi are restricted by the Second 
Law. For the sake of generality we do not assume any kind of reciprocity here. We 
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decouple the symmetric and antisymmetric parts introducing I = (J12 + Z21) /2 and 
k = (I12 — ^2i)/2. The entropy production is nonnegative if 

(38) h > 0, l 2 > and lj 2 - I 2 > 

Now the evolution equations (j3"6"|) - (l3"?l) can be written equivalently as 

(39) d t a = kB + hA + lB, 

(40) d t b = -kA + lA + i 2 B. 

3.3.1. Remark on dissipation potentials. We may introduce dissipation potentials 
for the dissipative part of the equations, if the condition of their existence is sat- 
isfied. Dissipation potentials are the children of variational principles, they are 
artificially added to a set of reversible equations of variational origin to generate 
some kind of dissipative effects. In our case there is no need of this assumption, 
our construction gives the most general dissipative system without any further ado. 
Moreover, here it is clear what belongs to the dissipative part and what belongs to 
the nondissipative part of the evolution equations. The terms with the symmetric 
conductivity contribute to the entropy production and the terms from the skew 
symmetric part do not. On the other hand, there is no need of potential construc- 
tion, as we are not looking for a variational formulation. Moreover, the symmetry 
relations are not sufficient for the existence of dissipation potentials in general, as 
we have emphasized previously. In case of constant coefficients (strict linearity), 
the dissipation potentials always exist for the dissipative (symmetric) part. 

3.3.2. Remark on the reciprocity relations. The reciprocity relations are the main 
results of the great idea of Lars Onsager connecting fluctuation theory to macro- 
scopic thermodynamics [601 ETJ [62, 63, 64J. As it was written by Onsager himself 
on the validity of his result: "The restriction was stated: on a kinetic model, the 
thermodynamic variables must be algebraic sums of (a large number of) molecular 
variables, and must be even functions of those molecular variables which are odd 
functions of time (like molecular velocities)" |63j . The Casimir reciprocity relations 
are based on microscopic fluctuations, too [64] . We do not have such a microscopic 
background for most of internal variables. E.g. in the case of damage the internal 
variables are reflecting a structural disorder on a mesoscopic scale. The relation 
between thermodynamic variables and the microscopic structure is hopelessly com- 
plicated. On the other hand, the Onsagerian reciprocity is based on time reversal 
properties of corresponding physical quantities either at the macro or at the micro 
level. Looking for the form of evolution equations without a microscopic model, 
we do not have any information on the time reversal properties of our physical 
quantities neither at the micro- nor at the macroscopic level. Therefore, we can 
conclude that lacking the conditions of the Onsagerian or Casimirian reciprocity 
gives no reasons to assume their validity in the internal variable theory. 

Let us observe the correspondence of evolution equations for internal variables 
with the reciprocity relations by means of a few simple examples. 

3.3.3. Example 1: Internal variables. Let us consider materials with diagonal con- 
ductivity matrix L (/ = 0, k = 0). It is clear that the Onsagerian reciprocity 
relations are satisfied, and we return to the classical situation with fully uncoupled 
internal variables: 

d t a — l\A, 
d t b = l 2 B. 

In this case the evolution equations for dual internal variables a and b are the same 
as in the case of single internal variable. 
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3.3.4. Example 2: Dynamic degrees of freedom. We now assume that all conduc- 
tivity coefficients are constant, and their values are h = I — 0, k — 1. This means 
that I12 = —I21, i.e., the Casimirian reciprocity relations are satisfied. 

For simplicity, we consider a specific decomposition of the entropy density into 
two parts, which depend on different internal variables 

s(a, dia, b, dib) — —K(b) — W(a, did). 

The negative signs are introduced taking into account the concavity of the en- 
tropy. Then the thermodynamic forces are represented as 

A = d a W -di{d dia W), B = d b K = K'(b), 

and Eqs. flM])-®]) are simplified to 

(41) d t a = B = K'(b) 

(42) d t b = -A + l 2 B = -d a W + d l (d dza W) + l 2 K'{b). 

One may recognize that the obtained system of equation corresponds exactly to a 
Hamiltonian form with the last term of ([4"2"|) as added dissipation. Here the entropy 
density s plays the role of a Hamiltonian density, a is a Lagrangian variable and b 
corresponds to the (slightly generalized) conjugated moment. The transformation 
into a Lagrangian form is trivial if K is quadratic K (b) = ^— , where m is a constant. 
Then the whole system corresponds to the general structure of dynamic degrees of 
freedom as one can compare to Eq. (5.14) in |52j with the Lagrangian 

£(<9 t a, a, dia) = m W(a,dia), and D(a, did) — ^^(<9 t a) 2 

as dissipation potential. Moreover, the entropy flux density (|34|) in case of our 
special conditions can be written as 

J i = -d dia sK'{b)+Zl 

and one can infer that natural boundary conditions of the variational principle 
related to the above Lagrangian correspond to the condition of vanishing entropy 
flow at the boundary, with 3o = 0*. 

Therefore, the variational structure of internal degrees of freedom is recovered 
in the pure thermodynamic framework. The thermodynamic structure resulted 
in several sign restrictions of the coefficients, and the form of the entropy flux is 
also recovered. The natural boundary conditions correspond to a vanishing extra 
entropy flux. 

3.3.5. Example 3: Diffusive internal variables [65, 66, 67 . Now we give an addi- 
tional example to see clearly the reduction of evolution equations of internal degrees 
of freedom to evolution equations for internal variables and the extension of the later 
to the previous one. 

We keep the values of conductivity coefficients (i.e., ii = / = 0, k = 1), but 
assume that both K and W are quadratic functions 

K(b) = ^b 2 , W(a,d ia ) = ^(d ia ) 2 , 

where a and (5 are positive constants according to the concavity requirement. 
In this case, the evolution equations (|4*T|) - (|4*2"|) reduce to 

(43) d t a = B = K'(b) = /3b, 

(44) d t b = ~d a W + d l (d 9ia W) + l 2 K'(b) = ad ll a + l 2 (3b. 
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Putting b from Eq. |g3J) into Eq. we have 

1 / 2 

(45) —r,dtta d t a = dua. 

ap a 

which is a Cattaneo-Vernotte type hyperbolic equation (telegraph equation) for the 
internal variable a. This can be considered as an extension of a diffusion equation 
by an inertial term or as and extension of a damped Newtonian equation (without 
forces) by a diffusion term. 

3.3.6. Discussion. In the framework of the thermodynamic theory with dual weakly 
nonlocal internal variables we are able to recover the evolution equations for internal 
degrees of freedom. 

We have seen that the form of evolution equations depend on the mutual interre- 
lations between the two internal variables. In the special case of internal degrees of 
freedom, the evolution of one variable is driven by the second one, and vice versa. 
This can be viewed as a duality between the two internal variables. In the case of 
pure internal variables of state, this duality is replaced by self-driven evolution for 
each internal variable. The general case includes all intermediate situations. 

It is generally accepted that internal variables are "measurable but not control- 
lable" (see e.g. [IE]). Controllability can be achieved by boundary conditions or 
fields directly acting on the physical quantities. We have seen how natural boundary 
conditions arise considering nonlocality of the interactions through weakly nonlocal 
constitutive state spaces. 

As we wanted to focus on generic inertial effects, our treatment is simplified from 
several points of view. Vectorial and tensorial internal variables were not considered 
and the couplings to traditional continuum fields result in degeneracies and more 
complicated situations than in our simple examples. 

It is important to remark, that skew symmetric couplings are not always re- 
lated directly to inertial effects and indicate two directions, one into mechanics and 
one into thermodynamics, where our method can be generalized. Let us mention 
here the related pioneering works of Verhas, where skew symmetric conductivity 
equations appear in different inspiring contexts |69[ 170] . 

Finally, let us mention that the idea of constructing a unified theoretical frame 
for reversible and irreversible dynamics has a long tradition. The corresponding 
research was not restricted to the case of internal variables and was looking for a 
classical Hamiltonian or a generalized variational principle that would be valid for 
both dissipative and nondissipative evolution equations (see, e.g., [71~ 1 \72 \ \73 \ 174] 
and the references therein). 

4. Classical Irreversible Thermodynamics 

In this section we investigate the evolution of extensive physical field quantities 
whose equation of motion is a balance and we will see that for first order weakly 
nonlocal state spaces the thermodynamic force-flux system is a consequence of the 
nonnegative entropy production. Treatment of Extended Irreversible Thermody- 
namics and further details can be found in [36] . 

A balance type evolution equation for a conserved quantity a can be given in a 
general form as 

(46) d t a+di? = 0. 

Here i 6 {1,2,3} denotes the spatial coordinates. The conserved quantity can 
be a Descartes product of densities of extensives of any tensorial order e.g. a = 
(p, e,p J , ...), where p is the mass density, e is the internal energy density, p 1 is the 
momentum density. Our continuum is at rest, j l are the corresponding fluxes, e.g. 
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'f = (pv l , pev 1 + jl, pv^v 1 + P JI , ...), where the first term is the current of the mass, 
the second is the convective and conductive current of internal energy and the third 
is the current of momentum. Therefore with this convenient notation we introduce 
only those free indices that are important from the point of view of the balance 
structure of the evolution equation. For the balances of particular real physical 
quantities there are peculiarities that we do not introduce here. For example for a 
continuum at rest traditionally there is no mass flux as a consequence of barycentric 
velocities, the source terms can play an important role (for internal energy, chemical 
production, etc.). Those peculiarities introduce further constraints and additional 
terms in the final equations, that we do not consider in this general calculation to 
show the core structure of classical irreversible thermodynamics. Moreover in this 
case we restrict ourselves for a continuum at rest. Some consequences of the motion 
of the continua is considered and investigated at the next section. We assume here 
a first order weakly nonlocal state space. 
Therefore in this case 

- the basic state space is spanned by a, 

- the constitutive state space is spanned by (a, Qja), 

- the constitutive functions are s, J 1 and j\ 

The above evolution equation is not completely arbitrary, but restricted by the 
Second Law of thermodynamics ([2]). Let us introduce a Lagrange-Farkas multiplier 
A for the evolution equation (|46[) and apply Liu procedure to determine the form 
of the constitutive functions, required by the entropy inequality: 

< <9 t s(a, d l a) + d l J\a, 9 4 a) - A (d t a + 9jf (a, 9* a)) 

= d a s d t a + d dtSL s d it a + d a J l dia + d dja J l d^a 

-A (d t a + dgfdia + d dja j l d tj aj 

= (d a s - Xj dta + d 9iSL s d it a + (d & .P - XdSj <%a 

(47) + (d dja J 1 - Xd d] J^ %a. 

The underlined partial derivatives of the basic field are not in the constitutive 
space, they are independent algebraic quantities and span the process direction 
space. The multipliers of those terms are zero according to Liu's theorem and give 
the Liu equations: 

(48) d t a : 9 a s = A, 

(49) d lt a : d 9lSL s = 0\ 

(50) %a : 8 d] J l =Xd d] J\ 

The first equation determines the Lagrange-Farkas multiplier and the second one 
shows that the entropy should be local, independent of the gradient of a. Therefore 
one can solve the third equation as 

(51) J l (a,d ia ) = d a s(a) +5 l (a), 

where we have denoted the variables of the constitutive functions and we have 
introduced the local residual entropy flux % l . Several authors suggest this kind of 
additive supplement to the classical entropy flux (e.g. the K vector of Muller [75]). 
This is a complete solution of the system (|48l) -(|50 |) . The dissipation inequality is 

(52) OKda'+fdiidJ). 

Assuming that the residual entropy flux is zero Z 1 = l , the entropy inequality 
reduces to the usual flux-force system of Classical Irreversible Thermodynamics, 
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where the thermodynamic forces are the gradients of the intensives and the ther- 
modynamic fluxes are identical to the the fluxes of the extensives from the balances. 

Flux: j\ Force: di(d & §). 

The classical solution of the above inequality is that the fluxes are proportional 
to the forces: 

(53) j* = t lk d k (d a s), 

where is symmetric and positive definite. Therefore balances of extensives are 
reduced to transport equations of the form: 

(54) Sta + d^L^cU)) =0. 

We have seen that in our calculations the classical form of the entropy production 
(the first term in the inequality above) and the classical form of the entropy flux 
were consequences of the Second Law, with the assumption of first order nonlocal- 
ity A rigorous treatment showed that the seemingly intuitive steps of irreversible 
thermodynamic modelling are well supported and explained by Liu procedure. We 
can see that the second part of the usual phenomcnological formulation of the lo- 
cal equilibrium hypothesis 72] is a consequence of the assumption of first order 
nonlocality: the state functions are those in equilibrium due to the local entropy. 



5. One component fluids - second order nonlocal in the density 

Up to know we have investigated evolution equations of internal variables and the 
general structure of Classical Irreversible Thermodynamics in a somewhat abstract 
manner. In this section we analyze a more particular and less abstract example 
where the physical meaning of the variables in the basic state space is well known. 
Our example is a one component heat conducting fluid where the constitutive state 
space for the energy and the velocity fields are first order, but second order in 
the density. We have seen that in case of first order weakly nonlocal state spaces 
our method gives the well known classical structure with a somewhat simplified 
manner, where the number of independent assumptions are reduced. E.g. the form 
of the entropy flux is calculated and not postulated, as we have demonstrated in 
the previous section. However, the second order nonlocality in the density variables 
leads to a surprising result and gives the viable family of Korteweg fluids, those 
that are compatible with the Second Law. One can find more details in [76l [77] 
and with alternate methods in (751 E] • 

5.1. Fluid mechanics in general. The basic state space of one-component fluid 
mechanics is spanned by the mass density p, the velocity v l and the energy density 
e of the fluid. Hydrodynamics is based on the balance of mass, energy and mo- 
mentum |71j . In classical fluid mechanics the constitutive space, the domain of the 
constitutive functions, is spanned by the basic state space (p, v l , e) and the gradient 
of the velocity diV 3 and the temperature diT. The pressure/stress tensor and the 
flux of the internal energy are the constitutive quantities in the theory. 
The balance of mass can be written as 

(55) dtp + =0. 

There is only a convective flux for mass. The balance of momentum, i.e. the 
Cauchy equation, is 

(56) d t {pv l ) + dj (p ij + /wV) = 0\ 

where pv l is the momentum density and P lJ is the pressure tensor, the conductive 
flux of the momentum. 
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The balance of energy is given as 



(57) 8 t e + dM" + zv l ) = 0, 

Here e is the total energy density and q" 1 is the conductive flux of the energy. The 
Second Law requires that the production of the entropy is nonnegative in insulated 
and source-free systems, therefore source terms are not considered in the balances. 
Hence in our case there is no mass production, the external forces are absent and 
the energy is conserved. We do not need the concept of the internal energy yet, 
first we analyse the consequences of the Second Law with the balance of the total 
energy. As we are dealing with a fluid, it is convenient to separate the conductive 
and convective entropy fluxes. Hence |(5J) will be written as 

(58) d t s + d t (J i + sv l ) > 0. 

Therefore in this case 

- the basic state space is spanned by (p, v l , e), 

- the constitutive state space is spanned by (p, dip, dijp, v l , diV 3 , e, die), 

- the constitutive functions are s 1 J l ,q l and P lJ . 

It is somewhat convenient to introduce the relative velocity v l as basic state 
variable (instead of the momentum) and the gradient of the energy as constitutive 
state variable (instead of the temperature gradient). As we have a second order 
weakly nonlocal extension in the mass density we need the gradient of (I55j) as a 
further constraint in the entropy inequality 

(59) d it p + d ij (pv j )=Q i , 

We introduce the Lagrange- Farkas multipliers A, A l ,Ti,j for the balances (f55l) . 
([55]). (|5SD and d57]) respectively. 

Now we apply Liu procedure, with the method of Lagrange-Farkas multipliers 
as in the previous sections 

< d t s + dj l + diisv*) - A (d t p + d^pv 1 )) - A* (d lt p + - 

(60) ti (d t {pv l ) + dj (P lJ + pwV)) - 7 (dte + di(qt + ev 1 )) . 

Developing the partial derivatives of the constitutive functions gives: 

dpsdtp + ddipsdup + do ijP sd U jP + d v isd t v l + d a . vi sd t jV % + d e sd t e+ 
d die sd ti e + d p J l dip + d djP J l d lj p + d ajkP J l d ljk p + d vj J l d. t v J + d dkV3 J l d jk v t + 

d e J l d l e + d dje J l d lJ e+ 
d t {sv l ) ~ A (d t p + pd t v l + v%p) - 
A 1 (pup + pd l] v 3 + v 3 dijp + dipdjv 3 + d^pdiV 3 ) - 

ti (pd t v l + v l d tP + pv'djv 3 + poldjV* + v^v'djp + il.l'-'i'.p + d 9kP P l3 d lkP + 
d dklP P l3 d jklP + d v uP ij d jV k + d aivk P ij d jl v k + d e P ij d je + d dke P ij d jk ej - 
7 (d t e + ediV 1 + v l d l e + d p q L d l p + i)„ ,,'/'('',,/' • 
ddj kP <?dijkP + d v jq l diV° + d dkV3 q l d lk v° + d e q l die + dg^fdjie) > 0. 
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A little rearrangement of the terms results in 

{d p s — A - f iv l )d t p + (d dzP s - k l )d tl p + d dijP sd U jp + [d v i§ - pf i )d t v i + 
ddivisdtjv 1 + (d e s - j)d t e + dd ie sd u e+ 
(d d]kP J l - ftdo^P 1 * - d d]hP it)d l]kP + 

(dg^r - tid dkVi P lj - id dkV ^d lk vi)d Jk v* - A l P d l y+ 

{d dj J l - ftd^pJ + *fd Bje ?)d ii e+ 
d p J l d l p + d d . p J l d io p + d v j Pd.v 3 + d e J i d i e+ 
d i (sv i )-X(pd i v i +v i d iP )- 
A l (pdijyi + v^dijp + dipdjv 3 + djpdiv^ — 
V, (/'/•'(), r ' + pr ■<),,■■ + r'v'i),p + dpP^djP + <h, „l"^) a p- 
d v kP ij djV k + <9 e P y d je ) - 
7 (ediV 1 + v i d i e + d p q l d t p + da^^djip + d vj q l d % v J + d e <f <%e) > 0. 

Then the Liu-equations follow as the multipliers of the members of the process 
direction space, the derivatives that are out of the constitutive space 



(61) 


d t p 


d p s — A — fV = 0, 






(62) 


dnP 


d diP s ~ A* = 0\ 






(63) 


dujP 


a, P s <)". 






(64) 


dtv 1 


d v i§ - pti = Qi, 






(65) 


dtjV* 


8g jV iS = Of, 






(66) 


d t e 


d e s - 7 = 0, 






(67) 


due 


d die S = 0\ 






(68) 


dijkP 


da k]P r = £id dkjP P li 


+ jd dkjP q\ 




(69) 


d ok v l 


d 6kVi p =f,d dkV >P l > 






(70) 


dije 


d dj J l ^ftfa^ + 


jd dje q\ 





As a consequence of (p5|) and |H7|) , the entropy density does not depend on 
dijp, djV 1 and c^e. ([CTjl. ((HJ), ([Ml) and give the Lagrange- Farkas multipliers in 
terms of the entropy derivatives. Therefore, from a thermodynamic point of view, 
the Lagrange-Farkas multipliers are the normal and generalized intensive variables 
[79]. Now, one can give a solution of ([6"5| - (|70[) as 

(71) P = d e sq l + £ (datpSBjV* + d d]P sd v l ) + -d vJ SP** + fr(p, d lP , v\ e), 

A p 

where the residual entropy flux Z l is an arbitrary function. Thus Liu's equations 
can be solved and yield the Lagrange-Farkas multipliers as well as restrictions for 
the entropy and the entropy flux. Applying these solutions of the Liu equations, 
the dissipation inequality can be simplified to the following form 

(72) < <7 S = dft + tfdiideb*)) + P ij di {d viS ) + 

djV j (^s + ed e s - pdps + (<9a, P s)^ + djV* f ( d d jP s) j . 
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Here we have introduced the specific entropy as s := s/p. It is worth to give this 
inequality, the entropy production in a more traditional form, without indices, too 



(73) V- J + q- V9 e (ps)- 



V<9 v s : P + 



P 2 

— (V • 5 Vp sI + V<9 Vp s) + (§ + ed e s - pd p s)l 



: Vv > 0. 



Here I denotes the second order unit tensor Sij, V is the gradient, V- is the 
divergence of the corresponding field quantity. To get the traditional form of the 
equation we should introduce the internal energy of the fluid as the difference of the 
total and kinetic energies and assume that the entropy function, does not depend 
on the total and kinetic energies independently, but only on the internal energy 
s = s(p, dip, e — pv 2 /2). Moreover, we want to define the entropy as an extensive 
quantity, therefore we require that the specific entropy depends on the specific 
internal energy e: 

/ 6 V 

(74) s{p,dip,v\e) = ps\p,-- —,dip 

From this form of the entropy function we get the Gibbs relation 

de = Tds + ^rdp - A l dd lP . 
P 2 

2 

where e = - — ^- is the specific internal energy. The temperature and pressure are 
defined by the customary partial derivatives of the entropy. The temperature can be 
connected both to the derivative by the total energy and the internal energy, because 
the corresponding derivatives are equal. A 1 is defined as the partial derivative of 
the entropy by the density gradient, like the traditional intensives. : 



de{ps) = ^, d vi s = ~, d p s = d dtP s = ^. 

With these quantities we can write the dissipation inequality and the entropy 
flux as 



(75) 0< a s =di^ + {q i -VjP ij )di-- 



1 

f 

1 ffyii (. , T P\ , a .s\ ri V 



- ^ -{p+ -j-d k (d dkp5 ) j 5) -j-dj (d dip5 ) ) dtf. 

(76) J 1 = {?-v j P* i )± + £(d BtP sd J vi+d BkP Sd k v i )+&. 

We can see that a flux of the internal energy is introduced as in case of the 
traditional, first order weakly nonlocal theories. However, the appearance of tem- 
perature inside the expression of the viscous pressure indicates the possibility of 
alternate, better definitions. 

Now we change the notation to a coordinate invariant one, introducing the nabla 
operator for the space derivatives as it is customary in fluid mechanics. In this case 
e.g. diV 1 = V • v and = Vv. In the pure mechanical, reversible case our 
thermodynamic force for mechanical interactions is zero. Therefore we introduce 
the nonlocal reversible pressure as 

(77) P r = ^ [(V • d Vp s - 2d p s) I + VdvpS] • 

If the pressure is equal to the reversible pressure, there is no dissipation, the 
theory is reversible (conservative). In case of a local entropy (independent of the 
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gradient of the density) then we obtain 

(78) Peu1c» := -T P 2 d pS (p)I, 

therefore, the corresponding equations are of the ideal Euler fluid, where p(p) = 
—Tp 2 d p s(p) is the scalar pressure function. Introducing the viscous pressure P v as 
usual, we can solve the dissipation inequality and give the corresponding Onsagerian 
conductivity equation as 

pv = p pr = Lqns . w 

Here 1*ons is a nonnegative constitutive function. Let us recognize that if s 
is independent of the gradient of the density, Lews is constant, and P v is an 
isotropic function of only Vv, then we obtain the traditional Navier-Stokes fluid 
(see e.g. [7T]). 

One can prove easily that the reversible part of the pressure is potentializable, 
i.e., there is a scalar valued function U such that 

(79) V • P r = pVU. 
U can be calculated from the entropy function as 

(80) U= V ■ (pdyps) - d p (ps). 

Therefore in case of reversible fluids the momentum balance can be written 
alternatively as 

(81) pv + V-P r = ^> v + VL> = 0. 

Let us give an interesting particular example of a weakly nonlocal fluid. 



5.2. Schrodinger-Madelung fluid. Here the entropy is defined as 



(82) WAVp,v).^g) 2 4 = ^(vm,f 2 , 

where v is a constant scalar. The corresponding reversible pressure is 



(83) 



P'- = -g(ApI + VV- 2V ^ ) 



where o denotes the tensorial/dyadic product, as mentioned before. The potential 
is 

(84) Us ^ = -T P { Af) --W- )=- HIT' 

where we introduced R — y/p to show that ([M]) is the quantum potential in the 
de Broglie-Bohm version of quantum mechanics (if v = h 2 /m 2 ) [501 181j . 
The entropy flux of the Schrodinger-Madelung fluid is 

(85) J SchM = -v • P r - J ( VpV ■ v + Vp • Vv) . 



5.2.1. Remark. The results of this subsection are strange with the traditional inter- 
pretation of objectivity and frame independence. That question will be discussed 
from a more general point of view in the next section. 
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5.2.2. Discussion. An important property of the Schrodinger-Madelung fluid is that 
if the motion of the fluid is vorticity free, V x v = 0, then the mass and momentum 
balances can be transformed into and united in the Schrodinger equation. Hence 
the balance of momentum (|56|) can be derived from a Bernoulli equation (in a given 
inertial reference frame). Defining a scalar valued phase (velocity potential) by 

v = ^VS, 
m 

we obtain the Bernoulli equation observing that the second part of (|8T|) is the 
gradient of 

ft dS V2 TT 

(86) ^m + Y- Us ^ = °- 

Then, introducing a single complex valued function ip :— Re lS that unifies R = 
yfp and S, it is easy to find that the sum of (|55| multiplied by ihe lS / (2R) and ([55]) 
multiplied by mRe form together the Schrodinger equation for free particles: 

<*>_ »S~5> 

It is remarkable that the structure of quantum mechanics appears in a classical 
thermodynamic approach without any explicit distinctive assumption related to the 
microscopic quantum world. In this sense the basic assumptions of our derivation 
are rather weak and very general. 



6. Summary and outlook 

In this paper we have shown that the Second Law of thermodynamics provides a 
general, uniform, rigorous and constructive method to investigate weakly nonlocal 
extensions of classical and nonclassical non-equilibrium thermodynamics. We have 
seen how one can generate evolution equations for internal variables, to understand 
the constitutive structure of Classical Irreversible Thermodynamics and to restrict 
the pressure tensor of ideal and viscous Korteweg fluids. 

In the suggested approach the choice of the constitutive state space is a phys- 
ical question. The number of the necessary constraints for a given constitutive 
state space was introduced intuitively, according to our calculational experience, 
always looking forward to an explicit solution of the Liu equations. Cimmelli in 
[82) introduces a different philosophy, requiring that the number of the constraints 
must equal the number of constitutive state variables. Additional gradients of the 
constraints are to be introduced in case of necessity. In general, the method given 
in [37] introduces less additional constraints, therefore the results are more easily 
applicable. On the other hand, the method proposed in |82) is more cumbersome 
and the interpretation of the consequences of the Second Law can be less straight- 
forward. However, it is based on a precise rule. 

One of the interesting observations was that we were able to derive some re- 
strictions to the reversible part of the evolution equations. There we have encoun- 
tered equations of Euler-Lagrange type in case of zero entropy production in the 
corresponding interaction (static Ginzburg-Landau, dynamic degrees of freedom, 
Bohm-potential for Korteweg pressure). In a sense we were able to derive Hamil- 
tonian variational principles, we have proved their existence from the Second Law 
of thermodynamics. 

There are several other classical and nonclassical weakly nonlocal equations of 
physics that could be investigated in this general frame. We have analyzed weakly 
nonlocal extensions of the heat conduction equations (both Fourier and Cattaneo- 
Vernotte) were researched, too. In this respect we have got that special kind of 
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internal variables, the so called current multipliers can represent some kind of non- 
local effects in the theory. They can lead to theoretical structures like the Guyer- 
Krumhansl equation of heat conduction [83] . This concept is originated in general 
extensions of the entropy flux jSU [85] . These possibilities have been investigated 
together with higher order weakly nonlocal extensions in case of generalized heat 
conduction and generalized Ginzburg-Landau equations in [861 150] . The flux hi- 
erarchy of extended irreversible thermodynamics and the origin of balance form 
evolution equations for internal variables was researched in |87[ 188] . 



On the other hand, we have researched phase separated multicomponent fluids 
and got a structure with natural instabilities (kind of nonlocal phase boundaries) 
similar to the Goodman-Cowin pressure [5"9l [§01 [9"T] . 

However, why should we restrict ourselves to weakly nonlocal extensions of the 
constitutive functions? Why do not we consider time derivatives in exploiting the 
Second Law? The answer to this question leads to one of the most fundamen- 
tal problems of physics, the question of objectivity, in particular, to the question 
of material frame indifference. In this paper we restricted ourselves to evolution 
equations and physical quantities given in an inertial reference frame. However, 
we know well, that the laws of physics are independent of an external observer. In 
particular, the constitutive functions can depend only on the material (e.g. on the 
motion of the continua) but not on the motion of an external observer. There are 
strong arguments that the usual formulation of the material frame indifference and 
the usual concept of objectivity is wrong 92 , 93, 94, 95j. These investigations show 
well, that the core of the problem is in the usual intuitive concept of the nonrcla- 
tivistic spacetime [55] . Therefore, we need more rigorous analysis of the kinematics 
of classical continua from this point of view and an absolute, frame independent 
formulation of the exploitation method of the Second Law. As regards the new 
foundations of finite deformation kinematics the investigations of Fiilop seem to 
settle a good starting point [97] . 

The results of those investigations are exploited in two key points in this paper. 
First of all, from the objectivity point of view, considering mathematical mod- 
els of nonrelativistic spacetime, the successes of gradient extensions can be well 
understood. The gradient of a physical quantity is a spacelike component of a 
four-covector (spacetime quantity) and therefore it is independent of an observer 
[98j : for inertial observers three- vectors are transforming by Galilei transforma- 
tion, but three-covectors are invariant. Therefore gradients are objective. On the 
other hand, for fluids we have started our investigations with an explicitly velocity 
dependent constitutive state space. That is excluded by the old concept of objec- 
tivity, because the relative three-velocities are frame dependent. However, it is not 
excluded according to the new concept of objectivity, because the four-velocities 
are independent of the frame. Our treatment of a one component fluid here is 
not a true objective treatment, from more than one point of view. However, we 
have done the first steps toward the development of objective exploitation methods 
of the Second Law, too. As the problem of objectivity is in a sense more easily 
understandable from the point of view of a relativistic spacetime model, we have 
investigated the nonequilibrium thermodynamics of special relativistic fluids [9"9"] . 
Moreover, we have started to investigate some possible practical consequences of 
incorporating our generalization of the objective time derivatives of rheology into 
a thermodynamic framework [1001 1101] . 

Finally let us emphasize how surprising is our result regarding fluid mechan- 
ics. The original observation of Madelung was that the Schrodinger equation can 
be transformed into an interesting fluid mechanical form [102] . The obervation of 
Bohm was that the same equation can be transformed into a Newtonian form [80] . 



20 



In both cases the Schrodinger equation is postulated, coming out of the blue. Here 
we have derived the fluid mechanical equivalent of the Schrodinger equation in a 
very general framework, from a surprisingly minimal set of assumptions, from the 
basic balances, the Second Law of thermodynamics and nonlocality of the inter- 
actions were used. The existence of such derivation is a disturbing fact that is 
hardly understandable from the traditional point of view regarding the relation of 
the irreversible macro- and reversible microworld. 
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7. Appendix - Farkas's lemma and some of its consequences 

Lemma 7.1. (Farkas) Let ^ be vectors in a finite dimensional vector space 
V, i = l...n, and S = {p G V*|p - a.; > 0, i = l...n}. The following statements are 
equivalent for all b G V: 

(i) p • b > 0, for all p G S. 

(ii) There are nonnegative real numbers X±, A„ such that b = Y]7—i A^a^. 
Proof: 

(ii) =^ (i) p • Y%=i x i a * = E"=i Aip • a, > if p G S. 

(i) => (ii) Let us consider a maximal, linearly independent subset ai, ...,a; of S. 
Let S = {y G V*|y • a; = 0,i = 1...1}. Clearly ^ So C S. 

If y G So then — y is also in So, therefore y • b > and — y ■ b > together. 
Therefore for all y G So it is true that y • b = 0. As a consequence b is in the 
linear subspace generated by {a^}, that is there are real numbers Ai, A/ such that 

b = Ei=l A * a i- 

Moreover, for all k G {1, ---J} there is a p& G V* such that p& • a& = 1 and 
Pfc-aj = if i 7^ k. Evidently, pfc G S for all k, therefore < Pfc-b = Pk-J2\=i = 
XiPk ■ SLi — Afc is valid for all k. Lastly, we can choose zero multipliers for the 
vectors that are not independent. ■ 

Remark 7.1. In the following the elements o/V* are called independent variables 
and V* itself is called the space of independent variables. The inequality in the first 
statement of the lemma is called objective inequality and the nonnegative numbers 
in the second statement are called Lagrange- Farkas multipliers. The inequalities 
determining S are the constraints. 

In the calculations an excellent reminder is to use Lagrange- Farkas multipliers 
similarly to the Lagrange multipliers in case of conditional extremum problems: 

n n 

p • b - ^ Ajp • aj = p • (b - ^ Aj • > 0, Vp G V* 

i=l i=l 

From this form we can read out the second statement of the lemma. 

Remark 7.2. The geometric interpretation of the theorem is important and graphic: 
if the vector b does not belong to the cone generated by the vectors a*, there exists 
a hyperplane separating b from the cone. 

7.1. Affine Farkas's lemma. This generalization of the previous lemma was first 
published simultaneously by A. Haar and J. Farkas as subsequent papers in the same 
journal, with different proofs [1031 1104] , Later it was rediscovered independently 
by others several times (e.g. |105l I106| ). 
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Theorem 7.2. (Affine Farkas) Let &i ^ be vectors in a finite dimensional vector 
space V and on real numbers, i = l...n and Sa = {p £ V*|p • > a i: i = l...n}. 
The following statements are equivalent for a b £ V and a real number (3: 

(i) p • b > [3, for all p £ Sa- 

(ii) There are nonnegative real numbers Ai, A„ such that b = EiLi an< ^ 
P < ELi A i^- 

Proof: 

(ii) => (i) p • b = p • X)" =1 A^aj = £" =1 A 4 p • a; > ^™ =1 ^ P- 

(i) =>■ (m) First we will show indirectly that the first condition of Farkas's lemma 

is a consequence of the first condition here, that is if (i) is true then p • b > 0, for 

all peS*. 

Thus let us assume the contrary, hence there is p' £ S, for which p' • b < 0. 
Take an arbitrary p £ Sa, then p + hp' £ Sa for all real nonnegative numbers k. 
But now (p + fcp') • b = p • b + kp' ■ b < [3, if k> That is a contradiction. 

Therefore, according to Farkas's lemma exist nonnegative Lagrange-Farkas mul- 
tipliers Ai, A„ such that b = Yh=i ^ a i- Hence (3 < inf peSA {P ' E"=i = 
inf P es A {J2iLi x iP ■ a = E"=i ^i a i- ■ 

Remark 7.3. The multiplier form is a good reminder again 

m mm 

(p-b-/?)-^A i (p-a i -ai) = p-(b-^A i -a i )-/? + ^A i a i > 0, VpeV*. 

i— 1 i— 1 i— 1 

Remark 7.4. The geometric interpretation is similar to the previous one, but with 
affine objects. If the line (one dimensional affine hyperplane) (h,(3) does not belong 
to the (affine) cone generated by (a^a^) there exists an affine hyperplane separating 
b from the cone. 

7.2. Liu's theorem. Here the constraints are equalities instead of inequalities, 
therefore the multipliers are not necessarily positive. 

Theorem 7.3. (Liu) Let a^ ^ be vectors in a finite dimensional vector space 
V and ai real numbers, i — l...n and Sl = {p £ V*|p • a, = on,i = l...n}. The 
following statements are equivalent for abeV and a real number (3: 

(i) P • b > (3, for all p £ S L , 

(ii) There are real numbers Ai, A„ such that 

n 

(88) b = ^A,a 2 , 

i=l 

and 

n 

(89) /3<^W 

i=l 

Proof: 

A straightforward consequence of the previous affine form of Farkas's lemma 
because Sl can be given in a form Sa with the vectors a^ and — a^, i = 1, ...,n: 
S L = {p £ V*|p ■ a, > cti,p ■ (-aj) > cti,i = l...n}. 

Therefore there are nonnegative real numbers A^, A^ and Aj~, A~ such, that 

b = Er=i(\ +a i - v a ») = Er=i(\ + - v) a * = e™=i A * a * and p < Er=i(\ +a i - 

X~ai). U 

Remark 7.5. The multiplier form is a help in the applications again 

n n n 

0< (p-b-/3)-^A i (p-a i -a i ) = p-(b-^A i -a i )-/? + ^A i a i , VpeV*. 

i— 1 i—l i—1 
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Remark 7.6. In the theorem with Lagrange multipliers for local conditional ex- 
tremum of differentiate function we apply exactly the above theorem of linear al- 
gebra after a linearization of the corresponding functions at the extremum point. 

Considering the requirements of the applications we generalize Liu's theorem to 
take into account vectorial constraints: 

Theorem 7.4. (vector Liu) Let A ^ in a tensor product V ® U of finite dimen- 
sional vector spaces V and U. Let a £ U and Sl = {p G V* |p • A = a}. The 
following statements are equivalent for a b £ V and a real number (3: 

(i) p • b > [3, for all p € Sl. 

(ii) There is a vector A in the dual o/U such that 

(90) b = A • A, 
and 

(91) P<X-a. 
Proof: 

Let us observe that we can get back the previous form of the theorem by in- 
troducing a linear bijection K : U — > R n , a coordinatization in U. Then by the 
vectors a,; := A • K* • e,, where e^ is the standard i-th unit vector in R™, and the 
numbers on :— K • a, we obtain that b = 53lLi = -A-S™=i ^i^* e i = A A where 
A := ELi ^ K * e *- ■ 

The previously excluded degenerate case of A = deserves a special attention. 
Then, of course, a = and 5^ = V*; this is, in fact, a degenerate case of all the 
previous theorems and evidently the following statements are equivalent for a b € V 
and a real number (3: 

(i) p-b>/3for allpe V*, 

(ii) b = and (3 < 0. 

In this case the assignment of the vector space of V* is based on the inequality (i). 
We encountered various degeneracies in the calculations of the previous sections. 

Remark 7.7. In continuum physics and thermodynamics the above algebraic theo- 
rems are applied to differential equations and inequalities. There the constraints are 
differential equations and therefore V is generated by derivatives of some constitu- 
tive functions in the differential equations. V* is spanned by the process directions, 
the derivatives of the fields in the constitutive state space that are not already there. 
The corresponding form of and \91\) are called Liu equation(s) and the dissi- 
pation inequality, respectively. 
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